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We present a systematic and selfconsistent analysis of four-quark charmonium states and applied 
it to study compact four-quark systems and meson-meson molecules. Our results are robust and 
should serve to clarify the situation of charmonium spectroscopy above the threshold production of 
charmed mesons. 
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Understanding of charmonium spectroscopy is chal- 
lenging for experimentalists and theorists alike. Charmo- 
nium has been used as the test bed to demonstrate the 
color Fermi-Breit structure of quark atoms obeying the 
same principles as ordinary atoms [l|. Its nonrelativis- 
tic character (v/c « 0.2 — 0.3) gave rise to an amazing 
agreement between experiment and simple quark poten- 
tial model predictions as cc states 0. The opening of 
charmed meson thresholds was expected to modify the 
trend in the construction of quark- ant iquark models. In 
the adiabatic approximation meson loops were absorbed 
into the static interquark potential. Thus, close to the 
threshold production of charmed mesons models required 
of an improved interaction Q. The corrections intro- 
duced to the quark-antiquark spectra explained some de- 
viations observed experimentally Q. 

Since 2003, with the discovery of several states in 
the open charm sector, we have witnessed a growth of 
puzzling new mesons, being D*j(2317), D s j (2460) and 
£>q(2308) the most prominent examples. Later, several 
new states have joined this exclusive group either in the 
open-charm sector, D s j(2860) and D s j(3040), or in the 
charmonium spectra, like the well established X(3872) 
and F(4260), Z(3930), X(3940), F(3940), X(4008), 
X(4160), X(4260), F(4350), and y(4660). In addition, 
the Belle Collaboration has reported the observation of 
similar states with non-zero electric charge: the Z(4430), 
the Zi(4040) and the Z 2 (4240) that have not yet been 
confirmed by other experiments and remain somewhat 
controversial Q. These new states do not fit, in general, 
the simple predictions of the quark-antiquark schemes 
and, moreover, they overpopulate the expected number 
of states in (simple) two-body theories. This situation 
is not uncommon in particle physics. For example, in 
the light scalar-isoscalar meson sector hadronic molecules 
seem to be needed to explain the experimental data @- 
8]. Also, the study of the NN system above the pion 
production threshold required new degrees of freedom to 
be incorporated in the theory, either aspions or as ex- 
cited states of the nucleon, i.e., the A 
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cussion suggests that charmonium spectroscopy could be 
rather simple below the threshold production of charmed 
mesons but much more complex above it. In particu- 
lar, the coupling to the closest (cc)(nn) system, referred 
to as unquenching the naive quark model [ll|, could be 
an important spectroscopic ingredient. Besides, hidden- 
charm four-quark states could explain the overpopulation 
of quark-antiquark theoretical states. Thus, the new ex- 
perimental discoveries are offering exciting new insights 
into the subtleties of the strong interaction. 

In an attempt to disentangle the role played by multi- 
quark configurations in the charmonium spectroscopy we 
have obtained an exact solution of the four-body prob- 
lem based on an infinite expansion of the four-quark wave 
function in terms of hyperspherical harmonics [l2T ] . The 
method is exact but is not completely adequate to study 
states that are close to, but below, the charmed meson 
production threshold. Such states are called molecular, 
in the sense that they can be exactly expanded in terms of 
a single singlet-singlet color vector. Close to a threshold, 
methods based on a series expansion fail to converge since 
arbitrary large number of terms are required to deter- 
mine the wave function. From our analysis, we concluded 
that those four-quark states with two different asymp- 
totic physical thresholds (as it is the case of the cenn 
system that may split either into a (cc)(nn) or (cn)(nc) 
two-meson states) can hardly present a bound state since 
the interaction between any pair of quarks contributes to 
the energy of one of the two physical thresholds. How- 
ever, we observed that the root mean square radius of 
a few channels did not grow in the same manner as in 
those channels clearly converging to an unbound two- 
meson threshold. Instead, their radius remained stable 
and their energy did not cease slightly decreasing. 

For this reason, we have used a different technique 
that we developed when studying baryon spectra with 
screened potentials and that showed to be very powerful 
close to a threshold [13[ . In this case, the hyperspherical 
harmonic expansion of the wave function was computa- 
tionally very expensive. Instead, we solved the Faddeev 
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equations for negative energies using the Fredholm de- 
terminant method that permitted us to obtain robust 
predictions even for zero-energy bound states. For the 
charmonium the situation is similar but simpler. Simi- 
lar because we are working on the region where meth- 
ods based on infinite expansions are inefficient, but sim- 
pler since it is a two-body problem, the scattering of two 
mesons. Thus, we solve the Lippmann-Schwinger equa- 
tion looking for attractive channels that may contain a 
meson-meson molecule. In order to account for all ba- 
sis states we allow for the coupling to charmonium-light 
two-meson systems. With this method we circumvent 
the uncertainties associated to the slow convergence of 
the hyperspherical harmonic method for large grand an- 
gular momenta. 

When we consider the system of two mesons Mi and 
M2 (Mi — D,D*) in a relative S— state interacting 
through a potential V that contains a tensor force then, 
in general, there is a coupling to the MiM 2 D— wave and 
the Lippmann-Schwinger equation of the system is 

poo 

tf"°" (p, p"; E) = Vf°" (p, p") + £ / P ' 2 dp' 

xV^' s '(p,p')- — ff w y,p";iS),(l) 

J E — p' /2/i + le J 

where t is the two-body amplitude, j, i, and E are the an- 
gular momentum, isospin and energy of the system, and 
is, £'s', £"s" are the initial, intermediate, and final orbital 
angular momentum and spin; p and \i are the relative 
momentum and reduced mass of the two-body system, 
respectively. In the case of a two D meson system that 
can couple to a charmonium-light two-meson state, for 
example when DD is coupled to J/^Slo, the Lippmann- 
Schwinger equation for DD scattering becomes 

poo 

Yj Yj J P^Pl^-jf^ (P<*iPl) 

y=D~D*,J/Hru ty=0,2 

x G 7 (£;p 7 ) t^ s f t ^( Pl ,p^E), (2) 

with a,/3 = D~D* , J/^uj. For bound-state E < that 
the singularity of the propagator is never reached, we can 
neglect ie in the denominator. By changing variables, 

where b is a scale parameter, and the same for p(p a ) and 
p"(pp). Replacing the integral from —1 to 1 by a Gauss- 
Legendre quadrature, we obtain a set of linear equations. 
If a bound state exists at an energy Eb, the matrix deter- 
minant is zero. We took the scale parameter b of Eqs. ^ 
and J2]) to be b = 3 fm _1 and used a Gauss-Legendre 
quadrature with TV = 20 points. 

We have consistently used the same interacting Hamil- 
tonian to study the two- and four-quark systems to guar- 



antee that thresholds and possible bound states are eigen- 
states of the same Hamiltonian. Such interaction con- 
tains a universal one-gluon exchange, confinement, and 
a chiral potential between light quarks [14]. We have 
solved the coupled channel problem of the DD, DD*, 
and D*D .In all cases we have included the coupling to 
the relevant (cc)(nn) channel (from now on denoted as 
J/^uj channels). 

As we study systems with well-defined C— parity, let 
us comment on the DD system. Since neither DD 
nor DD* are eigenstates of C— parity, it is necessary to 
construct the proper linear combinations. Taking into 
account that C(D) = D and C(D*) = —D (with a rela- 
tive minus sin between them), it can be found that [17] : 

D 1 = ^=[dD* +DD*) (4) 

and 

D 2 = ±=(dD* -DD*) (5) 

are the eigenstates corresponding to C = —1 and C — 
+1, respectively. This does not depend on the quan- 
tum numbers of the system because D and D* are not a 
particle-antiparticlc pair. Once the C— parity of a DD 
state is fixed, its isospin is also determined. In particu- 
lar, for a DD S— wave state, positive C— parity requires 
isospin 0, while negative C— parity implies isospin 1. 

Table Q] and Figs. Q] and [5] summarize our results. We 
have specified the quantum numbers and plotted the 
Fredholm determinant of the attractive channels. The 
rest, not shown on the table, are either repulsive or have 
zero probability to contain a bound state or a resonance. 

Let us remark that, of all possible channels, only a 
few are attractive. Of the systems made of a particle 
and its corresponding antiparticle, the J PC (I) = ++ (0) 
channel is always attractive. In general, the coupling 
to the r\ c r\ channel reduces the attraction, but there is 
still enough attraction to expect a resonance close and 
above the threshold. This channel is much more attrac- 
tive for the D*D system than for DD, thus, in the lat- 
ter one could expect a wider resonance. It is easy to 
explain the reason for such a close-to-bind situation with 
these quantum numbers. They can be reached from a 
two-meson system without explicit orbital angular mo- 
mentum, while through a simple cc pair it needs a unit 
of orbital angular momentum. Similar arguments were 
used to explain the proliferation of light scalar-isoscalar 
mesons The most attractive channel in the DD 

case is the J p (I) = l ++ (0) and can be explained as 
before, except the unity of intrinsic spin due to the D* 
meson. A simple calculation of the DD system (Eq. JT])) 
indicates that the J PC '{I) = l ++ (0) and 1+~(1) are 
degenerate. It is the coupling to the J/^w (Eq. ([2])) 
that breaks the degeneracy to make the l ++ (0) more at- 
tractive. The isospin 1 channel becomes repulsive due 
to the coupling to the lightest channel that includes a 
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FIG. 1: Fredholm determinant for the J PC (I) = l ++ (0) DD 
system. Solid (dashed) line: results with (without) coupling 
to the J/^lo channel. 
0.8 



TABLE I: Attractive channels for the two D— mesons system. 
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pion. Then, the existence of meson-meson molecules in 
the isospin one DD channels can be discarded. Using 
the coupling to the J/^oj, not present in the calcula- 
tions at the hadronic level of [HI, [l6[ , we obtain a bind- 
ing energy for the J PC (I) = l ++ (0) in the range 0—1 
MeV, in good agreement with the experimental measure- 
ments of X(3872) (see Fig. [TJ. This result supports the 
analysis of the Belle data on B — > K + J/^tt+it^ and 



FIG. 2: Fredholm determinant of the most attractive J PC (I) 
channels for the DD and D* D systems. 
0.8 
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B -> K + D 7fir° that favors the A" (3872) being a bound 
state whose mass is below the D°D threshold \V$. The 
existence of a bound state in the l ++ (0) DD channel 
would not show up in the DD system because of quantum 
number conservation. 

Finally, we have found that the J PC (I) = 2++ (0,1) 
D*D (see Fig. [5]) are also attractive due to the cou- 
pling to the J/^w and J/^p channels, respectively. This 
would give rise to new states around 4 GeV/c 2 and one 
experimental candidate could be the y(4008). In this 
case, such a resonance would also appear in the DD sys- 
tem for large relative orbital angular momentum, L = 2. 
A similar behavior can be observed in resonances pre- 
dicted for the A A system [H[. 

In all cases, being loosely bound states whose masses 
are close to the sum of their constituent meson masses, 
their decay and production properties must be quite dif- 
ferent from conventional qq mesons. Our calculation does 
not exclude a possible mixture of standard charmonium 
states in the channels where we have found attractive 
molecular systems. This admixture could explain some 
properties of the X(3872) [H,[20|. We would like to em- 
phasize the similarity of our results to those of Ref. [2l[ 
in spite of our different approach. Our treatment is gen- 
eral, dealing simultaneously with the two- and four-body 
problems and using an interaction containing gluon and 
quark exchanges instead of the simple two-body one-pion 
exchange potential of Ref. ;2l[. Nevertheless, we also 
concluded that the lighter meson-meson molecules are 
in the vector- vector and pseudoscalar-vector two-meson 
channels. Finally, let us remark that our approach could 
also be applied to the the cess sector. 

To summarize, we have performed the first system- 
atic analysis of four-quark hidden-charm states as com- 
pact states or meson-meson molecules. For the first 
time we have performed a consistent study of all quan- 
tum numbers within the same model. Our predictions 
robustly show that no deeply bound states can be ex- 
pected for this system. Only a few channels can be ex- 
pected to present observable resonances or slightly bound 
states. Among them, we have found that the DD sys- 
tem must show a bound state slightly below the threshold 
for charmed mesons production with quantum numbers 
J PC (I) = l ++ (0), that could correspond to the widely 
discussed A" (3872). Of the systems made of a particle 
and its corresponding antiparticle, DD and D* D , the 
J PC {I) = 0++(0) is attractive. It would be the only 
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candidate to accommodate a wide resonance for the DD 
system. For the D*D the attraction is stronger and 
structures may be observed close and above the charmed 
meson production threshold. Also, we have shown that 
the J PC {I) = 2++(0,l) D*7T channels are attractive 
due to the coupling to the J/vPw and J/^p channels. Due 
to heavy quark symmetry, replacing the charm quarks 
by bottom quarks decreases the kinetic energy without 
significantly changing the potential energy. In conse- 
quence, four-quark bottomonium mesons must also exist 
and have larger binding energies. An experimental effort 
in this direction will confirm or rule out the theoretical 
expectations. If the scenario presented here turns out to 
be correct, it will open a new interesting spectroscopic 
area. 

When this work was finished we learned that particular 
studies of some of the new charmonium states coincide 



with our theoretical predictions about the more attrac- 
tive quantum numbers [22j. 
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